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An argument is presented that the dark matter is axions, at least in part. It has three steps. 
First, axions behave differently from the other forms of cold dark matter because they form a 
rethermalizing Bose-Einstein condensate (BEC). Second, there is a tool to distinguish axion 
(~| , BEC from the other dark matter candidates on the basis of observation, namely the study of 

^ f the inner caustics of galactic halos. Third, the observational evidence for caustic rings of dark 

matter is consistent in every aspect with axion BEC, but not with the other proposed forms 
of dark matter. 
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One of the outstanding problems in science today is the identity of the dark matter of the 
universe . The existence of dark matter is imphed by a large number of observations, including 
^> \ the dynamics of galaxy clusters, the rotation curves of individual galaxies, the abundances of 

■ light elements, gravitational lensing, and the anisotropies of the cosmic microwave background 
• radiation. The energy density fraction of the universe in dark matter is 23%. The dark matter 
I must be non-baryonic, cold and collisionless. Particles with the required properties are referred 

■ to as 'cold dark matter' (CDM). The leading CDM candidates are weakly interacting massive 
particles (WIMPs) with mass in the 100 GeV range, axions with mass in the 10"^ eV range, 
and sterile neutrinos with mass in the keV range. 

Based on work by my collaborators and I over the past three or four years, I have come to 
, , the conclusion that the dark matter is axions, at least in part. The starting point is that cold 

. dark matter axions thermalize as a result of their gravitational self-interactions'^. When they 

I thermalize, they form a Bose-Einstein condensate. It may seem surprising that axions thermalize 

as a result of their gravitational self-interactions since gravitational interactions among particles 
are usually thought to be negligible. Gravitational interactions among particles are in fact almost 
always negligible but cold dark matter axions are an exception because the axions occupy in 
huge numbers a small number of states (the typical quantum state occupation number is 10^^) 
and those states have enormous correlation lengths (of order parsec to Gpc, today). 

Let us call F = l/r the axion thermalization rate. On time scales short compared to r, cold 
dark matter axions form a degenerate Bose gas described by a classical field equation. Their 
behavior is then indistinguishable from that of ordinary CDM except on length scales that are 
too short (10^^ cm or so) to be of observational interest. On times scales large compared to r, 
cold dark matter axions thermalize. The thermalization of a degenerate Bose gas is a quantum- 
mechanical entropy generating process, not described by classical field equations. On time scales 
larger than r the axion state, i.e. the state that most axions are in, tracks the lowest energy 
state available to them. The behaviour of such a rethermalizing axion BEC is different from 
that of ordinary CDM and the differences are observable. 



The thermaUzation of cold dark matter axions is discussed in detail in ref. It is found 
there that rethermalization of the axion BEC by gravitational self-interactions is sufficiently fast 
that the axions that are about to fall into a galactic potential well almost all go to the lowest 
energy state consistent with the angular momentum they have acquired from tidal torquing. 
That state is one of net overall rotation, implying V x v ^ where v{x, t) is the velocity field 
of the infalling dark matter. In contrast, ordinary cold dark matter (e.g. WIMPs and sterile 
neutrinos) falls in with an irrotational velocity field, V x -u = 0. The inner caustics of galactic 
halos are different in the two cases. If the dark matter falls in with net overall rotation, the 
inner caustics are rings whose cross-section is a section of the elliptic umbilic (-D-4) catastrophe, 
called caustic rings for short'^. If the velocity field of the infalling particles is irrotational, the 
inner caustics have a 'tent-like' structure which is described in detail in ref. ^ and which is quite 
distinct from that of caustic rings. Evidence was found for caustic rings of dark matter. The 
evidence is summarized in ref. It is shown in ref. ^ that the evidence for caustic rings of dark 
matter is precisely and in all respects consistent with the predictions of a rethermalizing axion 
BEC. 

The above is the gist of the argument. It is elaborated further in the three sections below. A 
few comments may be in order. One question is: what fraction of the dark matter must be axions 
to justify the evidence for caustic rings. We hope to comment on this soon. Another question 
is: to what extent does the evidence for caustic rings require the dark matter to be QCD axions 
as opposed to some other kind of axion-like particle(s). The evidence requires that a sizable 
fraction of the dark matter be identical bosons, whose number is conserved on cosmological time 
scales, and which are sufficiently cold and thermalize sufficiently fast that they form a BEC. 
Furthermore the BEC must rethermalize sufficiently fast that the particles go to a state of net 
overall rotation as they are about to fall into galactic potential wells. It happens that the QCD 
axion with mass of order 10~^ eV has all these properties and since it solves in addition the 
strong CP problem of the Standard Model of elementary particles, it is reasonable to assume 
that the dark matter is in fact QCD axions. However, there are many axion-like particles^ 
that can equally well reproduce the evidence for caustics rings. Furthermore, whether or not 
the particle in question is the QCD axion, the prediction of Bose-Einstein condensation and 
subsequent caustic ring formation is rather insensitive to the particle mass and therefore does 
not provide a good guide to it. The axion is being searched for as a constituent of the Milky 
Way haloM 

cis db particle radiated by the Sun^^. and in experiments that convert photons to 
axions and axions back to photons behind a wall^^^. 

Finally, many authors have proposed^^^IElti^at the dark matter is a Bose-Einstein condensate 

of particles with mass of order 10~^^ eV or less. When the mass is that small, the dark matter 

BEC behaves differently from CDM on scales of observational interest as a result of the tendency 

of the B EC to delocalize. Due to the Heisenberg uncertainty principle, a BEC has Jeans' length 
I15I2I16I 

= il6.Gpmr'^ = 1.02 . 10- cm ( '-^) ^ ( l^llll^) ^ , (1) 



where p is the BEC density and m the constituent particle mass. As mentioned earlier, this 
length scale is unobservably small in the QCD axion case. In contrast, when m ~ 10^^^ eV, the 
Jeans' length is of order 3 kpc and has implications for observation. It leads to a suppression 
of the dark matter density near the galactic center. This may be a remedy for the excessive 
concentration of dark matter near galactic centers seen in numerical simulatons of structure 
formation with ordinary cdmEI. 



1 Bose-Einstein condensation of cold dark matter axions 



Shortly after the Standard Model of elementary particles was established, the axion was postu- 
lated'^ to explain why the strong interactions conserve the discrete symmetries P and CP. For 
our purposes the action density for the axion field may be taken to be 

>Ca = -\d,^d^^ - \m\^ + A/ _ ... (2) 



m 



where m is the axion mass. The self-coupling strength is' 

A = ^/^i±43-0.35^ (3) 

in terms of the axion decay constant / and the masses rriu and of the up and down quarks. 
In Eq. (j2]), the dots represent higher order axion self- interactions and interactions of the axion 
with other particles. All axion couplings and the axion mass 

^^6-10-6 eV^^^^^ (4) 

are inversely proportional to /. / was first thought to be of order the electroweak scale, but its 
value is in fact arbitrarj^^. However, the combined limits from unsuccessful searches in particle 
and nuclear physics experiments and from stellar evolution require / ^ 3 • 10^ GeV^^. 

Furthermore, an upper limit / < 10^^ GeV is provided by cosmology because light axions 
are abundantly produced during the QCD phase transition^SHI. jj^ spite of their very small mass, 
these axions are a form of cold dark matter. Indeed, their average momentum at the QCD epoch 
is not of order the temperature (GeV) but of order the Hubble expansion rate (3 • lO'^ eV) then. 
In case inflation occurs after the Peccei-Quinn phase transition their average momentum is 
even smaller because the axion field gets homogenized during inflation. In addition to the cold 
axion population, there is a thermal axion population with average momentum of order the 
temperature. 

The non-perturbative QCD effects that give the axion its mass turn on at a temperature of 
order 1 GeV. The critical time, defined by m{ti)ti = 1, is ti ~ 2 • 10"^ sec (//lO^^ GeV)i Cold 
axions are the quanta of oscillation of the axion field that result from the turn on of the axion 
mass. They have number density 

^ llo^G^ j [-^J 

where a{t) is the cosmological scale factor. Because the axion momenta are of order ^ at time 
ti and vary with time as a{t)~^, the velocity dispersion of cold axions is 

~ 4 ^ (6) 

mti a{t) 

if each axion remains in whatever state it is in, i.e. if axion interactions are negligible. The 
average state occupation number of cold axions is then 

That Af is much larger than one tells us that the effective temperature of cold axions is much 
smaller than the critical temperature for Bose-Einstein condensation. 



Bose-Einstein condensation may be briefly described as follows: if identical bosonic particles 
are highly condensed in phase space, if their total number is conserved and if they thermalize, 
most of them go to the lowest energy available state. The condensing particles do so because, by 
yielding their energy to the remaining non-condensed particles, the total entropy is increased. 
Eq. ([7]) tells us that the first condition is overwhelmingly satisfied. The second condition is also 
satisfied because all axion number violating processes, such as their decay to two photons, occur 
on time scales vastly longer than the age of the universe. The only condition for axion BEC 
that is not manifestly satisfied is thermal equilibrium. 

Axions are in thermal equilibrium if their relaxation rate T is large compared to the Hubble 
expansion rate H{t) = ^. The relaxation rate F is given in the particle kinetic regime by 

Fr^na 6vM (8) 

where a is the relevant scattering cross-section. The particle kinetic regime is defined by the 
condition E << 6E where 6E is the energy dispersion of the particles. Cold dark matter axions 
are however mostly in the opposite regime: E >> 5E which we call the condensed regime. 
Thermalization in the condensed regime is discussed in detail in our recent paper We find 
that the relaxation rate of cold axions through their Xcp"^ self-interactions is of order^^Sl 

Ea ~ n . (9) 

Tx{t)/H{t) is of order one at time ti but decreases as t a{t)~^ oc a{t)~^ afterwards, implying that 
cold axions briefiy thermalize as a result of their A^^ interactions when they are first produced 
during the QCD phase transition but, after this brief period of thermalization, the axions are 
decoupled again. 

However the axions rethermalize later as a result of t heir gravitational self- interactions. 
Their relaxation rate by gravitational interactions is of order^^EI 

Eg ~ AttG n (10) 



where / ~ (jnbv) ^ is their correlation length. Tg{t)/H(t) is of order 5 ■ 10 GeV)3 
at time ti but grows as ta~^{t) oc a{t). Thus gravitational interactions cause the axions to 
thermalize and form a BEC when the photon temperature is of order 500 eV (//lO^^ GeV)2. 
Bose-Einstein condensation causes the axion correlation length to grow until it becomes of order 
the horizon. The growth in the correlation length causes the thermalization to accelerate; see 
Eq. (jlOp . When I is some fraction of t, Tg{t)/H{t) oc a(t)~'^t^, implying that thermalization 
occurs on ever shorter time scales compared to the Hubble time. The question is now whether 
axion BEC has implications for observation. 



2 Dark matter caustics 

The study of the inner caustics of galactic halos^^^^ provides a useful tool. An isolated galaxy 
like our own accretes the dark matter particles surrounding it. Cold collisionless particles falling 
in and out of a gravitational potential well necessarily form an inner caustic, i.e. a surface of 
high density, which may be thought of as the envelope of the particle trajectories near their 
closest approach to the center. The density diverges at caustics in the limit where the velocity 
dispersion of the dark matter particles vanishes. Because the accreted dark matter falls in and 
out of the galactic gravitational potential well many times, there is a set of inner caustics. In 
addition, there is a set of outer caustics, one for each outfiow as it reaches its maximum radius 
before falling back in. We exploit the catastrophe structure and spatial distribution of the inner 
caustics of isolated disk galaxies. 



The catastrophe structure of the inner caustics depends mainly on the angular momentum 
distribution of the infalling particles ^. There are two contrasting cases to consider. In the first 
case, the angular momentum distribution is characterized by 'net overall rotation'; in the second 
case, by irrotational flow. The archetypical example of net overall rotation is instantaneous rigid 
rotation on the turnaround sphere. The turnaround sphere is defined as the locus of particles 
which have zero radial velocity with respect to the galactic center for the first time, their outward 
Hubble flow having just been arrested by the gravitational pull of the galaxy. The present 
turnaround radius of the Milky Way is of order 2 Mpc. Net overall rotation implies that the 
velocity field has a curl, V x i; ^ 0. The corresponding inner caus tic is a closed tube whose 



cross-section is a section of the elliptic umbilic (-D_4) catastrophe'^I^. We call it a 'caustic ring', 
or 'tricusp ring' in reference to its shape. In the case of irrotational flow, V x iT = 0, the inner 
caustic has a tent-like structure quite distinct from a caustic ring. Both types of inner caustic 
are described in detail in refP. 

If a galactic halo has net overall rotation and its time evolution is self-similar, the radii of 
its caustic rings are predicted in terms of a single parameter, called jmax- Self-similarity means 
that the entire phase space structure of the halo is time independent except for a rescaling of all 
distances by R{t), all velocities by R{t)/t and all densities by I2U 22 | 23 | 6 | Yqj- definiteness, 
R{t) is taken to be the turnaround radius at time t. If the initial overdensity around which the 
halo forms has a power law profile^^I] 

''''^(4.)' . (11) 



Mi \Mi, 

where Mi and 6Mi are respectively the mass and excess mass within an initial radius rj, then its 
subsequent evolution is self-similar with R(t) cx fs+w. In an average sense, e is related to the 
slope of the evolved power spectrum of density perturbations on galaxy scales. The observed 
power spectrum implies that e is in the range 0.25 to O.SS'^Sl. The prediction for the caustic ring 
radii is (n = 1, 2, 3, .. )SE1 



40 kpc / i;i.ot 



n V220 km/s/ V0.18 



(12) 



where Vj-ot is the galactic rotation velocity. Eq.( [T2]) is for e = 0.3. The a„ have a small e 
dependence. However, the oc 1/n approximate behavior holds for all e in the range 0.25 and 
0.35. 

Observational evidence for caustic rings of dark matter with the radii predicted by Eq. (I12p 
was found in: the statistical distribution of bumps in a set of 32 extended and well-measured 
galactic rotation curves, the distribution of bumps in the rotation curve of the Milky Way, the 
appearance of a triangular feature in the IRAS map of the Milky Way in the precise direction 
tangent to the nearest caustic ring, and the existence of a ring of stars at the location of the 
second (n = 2) caustic ring in the Milky Way. The observational evidence for caustic rings of 
dark matter is summarized in ref. The recent improved measurement^^!] gf the rotation curve 
of our nearest large neighbor, the Andromeda galaxy, provides new evidence. The new rotation 
curve shows three prominent bumps at radii 10 kpc, 15 kpc and 29 kpc, whose ratios accord 
with Eq. (fT2]l . 

To reproduce the evidence for caustic rings of dark matter, the specific angular momentum 
distribution on the turnaround sphere should be given by 

R(f\^ 

£{h, t) = jmax fix (zxh) —— (13) 



where h is the unit vector pointing to a position on the turnaround sphere, z is the axis of rotation 
and jmax is the parameter that appears in Eq. (fT2]) . Eq. (fT3|) states that the turnaround sphere 



at time t rotates with angular velocity u) = ^^^z. Each property of the angular momentum 
distribution (|13|) maps onto an observable property of the inner caustics: net overall rotation 
causes the inner caustics to be rings, the value of jmax determines their overall size, and the 
time dependence given in Eq. ()13p is responsible for a„ oc 1/n. We now show that each of these 
three properties follows from the assumption that the infalling dark matter is a rethermalizing 
axion BEC. 



3 Three successes 



3. 1 Magnitude of angular momentum 

We make the standard assumption that the angular momentum of a galaxy is due to the tidal 
torque applied to it by nearby protogalaxies early on when density perturbations are still small 
and protogalaxies close to one another^^Sl^ 'pj^g amount of angular momentum galaxies typically 
acquire by tidal torquing can be reliably estimated by numerical simulation because it does 
not depend on any small feature of the mass configuration, so that the resolution of present 
simulations is not an issue in this case. The dimensionless angular momentum parameter 

, (14) 

GM-2 

where G is Newton's gravitational constant, L is the angular momentum of the galaxy, M its 
mass and E its net mechanical (kinetic plus gravitational potential) energy, was found 
have median value 0.05. In the caustic ring model the magnitude of angular momentum is given 
by jmax- The evidence for caustic rings implies that the jmax-distribution is peaked at jmax — 
0.18. The relationship between jmax and A is^ 



X - / 6 8 1 . 

^-y^3^ - • 

For jmax =0.18 and e in the range 0.25 to 0.35, Eq. ([15]) implies A =0.051. The excellent 
agreement between jmax and A gives further credence to the caustic ring model. Indeed if the 
evidence for caustic rings were incorrectly interpreted, there would be no reason for it to produce 
a value of jmax consistent with A. 



3.2 Net overall rotation 

Next we ask whether net overall rotation is an expected outcome of tidal torquing. The answer 
is clearly no if the dark matter is ordinary CDM. Indeed, the velocity field of ordinary CDM 
satisfies 

dv dv — » -* 

-^ir,t) = —if,t) + {vif,t) ■ V)vif,t) = -Vcj){f,t) (16) 

where (j){f, t) is the gravitational potential. The initial velocity field is irrotational because the 
expansion of the universe caused all rotational modes to decay away. Furthermore, it is easy to 
show^ that if V X {r = initially, then Eq. (jl6p implies V x w = at all later times. Since net 
overall rotation requires V x u 7^ 0, it is inconsistent with ordinary CDM, such as WIMPs or 
sterile neutrinos. If WIMPs or sterile neutrinos are the dark matter, the evidence for caustic 
rings, including the agreement between jmax 

and A obtained above, is fortuitous. 
Axions do not obey Eq. (jl6p because they form a rethermalizing By rethermalizing 

we mean that the thermalization rate is larger than the Hubble rate so that the axion state 
tracks the lowest energy available state. The compressional (scalar) modes of the axion field are 
unstable and grow as for ordinary CDM, except on length scales too small to be of observational 



interest"^. Unlike ordinary CDM, however, the rotational (vector) modes of the axion field 
exchange angular momentum by gravitational interaction. Most axions condense into the state 
of lowest energy consistent with the total angular momentum, say L = Lz, acquired by tidal 
torquing at a given time. To find this state we may use the WKB approximation because the 
angular momentum quantum numbers are very large, of order 10^*^ for a typical galaxy. The 
WKB approximation maps the axion wavefunction onto a flow of classical particles with the same 
energy and momentum densities. It is easy to show that for given total angular momentum the 
lowest energy is achieved when the angular motion is rigid rotation. Rigid rotation is therefore 
a prediction of tidal torque theory if the dark matter is axions. 

3.3 Self- similarity 

Because the axion BEC rethermalization rate is larger then the Hubble rate, most axions go 
to the lowest energy state consistent with the total angular momentum acquired from nearby 
inhomogeneities by tidal torquing. The time dependence of the specific angular momentum 
distribution on the turnaround sphere is then predicted. Is it consistent with Eq. (|13p ? In 
particular, is the axis of rotation constant in time? 

Consider a comoving sphere of radius S(t) = Sa{t) centered on the protogalaxy. As before, 
a{t) is the cosmological scale factor. S is taken to be of order but smaller than half the distance 
to the nearest protogalaxy of comparable size, say one third of that distance. The total torque 
applied to the volume V of the sphere is 

T{t)= d-^r 6p{f,t) fx {-V4>{f,t)) (17) 
Jv(t) 

where 5p{f,t) = p{r,t) — po{t) is the density perturbation. po{t) is the unperturbed density. 
In the linear regime of evolution of density perturbations, the gravitational potential does not 
depend on time when expressed in terms of comoving coordinates, i.e. (f){r = a{t)x,t) = <j){x). 
Moreover 6{f,t) = ^^^^(f-^ has the form 5{r = a{t)x,t) = a(t)6{x). Hence 

f(t) = po{t)a{t)^ f d?x 5{x) X X {-W^<i){x)) . (18) 
Jv 

Eq. (jlSp shows that the direction of the torque is time independent. Hence the rotation axis 
is time independent, as in the caustic ring model. Furthermore, since pQ (t) cx a{t)-^, T{t) 

2 5 2 2 

a{t) cx t3 and hence £{t) oc L{t) cx ts. Since R{t) a is+w, tidal torque theory predicts the 
time dependence of Eq. (fT3]) provided e = 0.33. This value of e is in the range, 0.25 < e < 0.35, 
predicted by the evolved spectrum of density perturbations and supported by the evidence for 
caustic rings. So the time dependence of the angular momentum distribution on the turnaround 
sphere is also consistent with the caustic ring model. 

In conclusion, the phase space structure of galactic halos predicted by tidal torque theory, if 
the dark matter is axions, is precisely and in all respects that of the caustic ring model proposed 
earlier on the basis of observations. The other dark matter candidates predict a different, more 
chaotic phase space structure for galactic halos. Although the QCD axion is best motivated, a 
broader class of axion-like particles behaves in the manner described here. 
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